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We propose novel mixed states in two qubits, "maximally 
entangled mixed states", which have a property that the 
amount of entanglement of these states cannot be increased 
further by applying any unitary operations. The property 
is proven when the rank of the states is less than 4, and 
confirmed numerically in the other general cases. The cor- 
responding entanglement of formation specified by its eigen- 
values gives an upper bound of that for density matrices with 
same eigenvalues. 



Entanglement (or inseparability) is one of the most 
striking features of quantum mechanics and an impor- 
tant resource for most applications of quantum informa- 
tion. In quantum computers, the quantum information 
stored in quantum bits (qubits) is processed by operat- 
ing quantum gates. Multi-bit quantum gates, such as the 
controllcd-NOT gate, are particularly important, since 
these gates can create entanglement between qubits. 

In recent years, it has attracted much attention to 
quantify the amount of entanglement, and a number of 
measures, such as the entanglement of formation Q and 
negativity ||,[|, have been proposed. When the system 
of the qubits is in a pure state, the amount of entangle- 
ment can be changed through the gate operations from 
zero of separable states to unity of maximally entangled 
states. The most quantum algorithms are designed for 
such ideal pure states. When the system is maximally 
mixed, however, we cannot receive any benefit of entan- 
glement in the quantum computation, since the density 
matrix of the system (unit matrix) is invariantly separa- 
ble under any unitary transformations or gate operations. 
Recently, a question about NMR quantum computation 
has been proposed 0, since the states in the vicinity of 
the maximally mixed state are also always separable as 
is the case of the present NMR experiments. 

In all realistic systems, the mixture of the density ma- 
trix describing the qubits is inevitably increased by the 
coupling between the qubits and its surrounding environ- 
ment. Therefore, it is extremely important to understand 
the nature of entanglement for general mixed states be- 
tween two extremes of pure states and a maximally mixed 
state. 

In this paper, we try to answer a simple question how 
much the increase of the mixture limits the amount of 
entanglement to be generated by the gate operation or 



unitary transformation. We propose a class of mixed 
states in bipartite 2x2 systems (two qubits) . The states 
in this class, maximally entangled mixed states, show a 
property of having a maximum amount of entanglement 
in a sense that the entanglement of formation (and even 
negativity) of these states cannot be increased further by 
applying any (local and/or nonlocal) unitary transfor- 
mations. The property is rigorously proven in the case 
that the rank of the states is less than 4, and confirmed 
numerically in the case of rank 4. The corresponding en- 
tanglement of formation specified by its eigenvalues gives 
an upper bound of that for density matrices with same 
eigenvalues. 

The entanglement of formation (EOF) 0] for a pure 
state is defined as the von Neumann entropy of the re- 
duced density matrix. The EOF of a mixed state is de- 
fined as Ep(p) = min^ i PiEp('ip) 1 where the minimum 
is taken over all possible decompositions of p into pure 
states p — X^-PilV'iXV'il- The analytical form for EOF in 
2x2 systems is given by || 



E , (p)= „ ( i±4Z^), 



id 



with H{x) being Shannon's entropy function. The con- 
currence C is given by 



C = max{0, Ai — A2 — A3 — A4}, 



(2) 



where A's are the square root of eigenvalues of pp in de- 
creasing order. The spin-flipped density matrix p is de- 
fined as 



(3) 



where * denotes the complex conjugate in the computa- 
tional basis. Since Ef is a monotonic function of C, the 
maximum of C corresponds to the maximum of Ef- 

The states we propose are that obtained by applying 
any local unitary transformations to 



M = pi|*-)(*-| +p 2 |00)(00| 

+ P3|* + )(* + |+P4|ll)(ll|, 



(4) 



where |* ± ) = (|01)±(10))/V2 are Bell states, and |00) and 
1 11) are product states orthogonal to \$ ). Here, p^s are 
eigenvalues of M in decreasing order (p\ > P2 > P3 > P4), 
and pi +P2+P3+P4 — 1- These include states such as 

p = Pl \<p-)(<p-\ + P2 |01)(01| 
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P3|<f + }(<I> + |+P4|10)<10|, 



(5) 



where 1$*) = (|00) ± |11))/V2 are also Bell states, and 
include such that obtained by exchanging | V E'~) <-> |\& + ), 
|00) <-> |11) in Eq. (g), or |$") <-> |$+), |01) ^ |10) 
in Eq. ([5]). Since entanglement is preserved by local 
unitary transformations, all these states have the same 
concurrence of 



C* = max{0, C*( Pi )} 

C*(pi) =P1-P3~ 2\/P2P4 



(6) 



The concurrence C* is maximum among the density ma- 
trices with the same eigenvalues, at least, when the den- 
sity matrices have a rank less than 4 (p^ = 0) . The proof 
is as follows: 

(1) Rank 1 case (p2 =P3 =Pa = 0). In this case, Eq. 
(Q) is reduced to M = \^~){^~\ which obviously has the 
maximum concurrence of unity. 

(2) Rank 2 case (p3=J>4 = 0). Any density matrices of 
two qubits (not necessarily rank 2) can be expressed as 



p = q\ip)(ip\ + (l-q) Pl 



scpj 



(7) 



where \ip) is an entangled state and p sep is a separable 
density matrix. Convexity of the concurrence [fffl implies 
that 

C{p) < qCm(i>\) + (1 - q)C( Pscp ) = qC{\i>m). (8) 

Since p sep is a positive operator, q is equal to or less than 
the maximum eigenvalue of p, and thus, 



C{p)<px. 



(9) 



The equality is satisfied when \ip) is a maximally entan- 
gled pure state and it is an eigenvector of p with the 
eigenvalue of p\ . The upper bound in Eq. (g) coincides 
with C* for p3=pi = 0. 

(3) Rank 3 case (p4 = 0). Any rank 3 density matrices 
can be decomposed into two density matrices as 



p = (1 - 3p 3 )p2 + 3p 3 P3, 

where eigenvalues of p2 are 

i l-3p3'l-3p 3 ' U,Uh 



(10) 



(11) 



and eigenvalues of p3 are {1/3, 1/3, 1/3, 0}. According to 
Lemma 3 in Ref. fl , 



Trp 2 < - p is separable. 
3 



(12) 



Since the purity of p% is 1/3, p% is always separable. 
Therefore, convexity of the concurrence implies that 
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C(p) < (l-3p 3 )C(p 2 ) <pi- P 3 - 



(13) 



FIG. 1. Numerically obtained maximum concurrence for 
random density matrices as a function of the participation 
ratio and C* (inset). 



Here, we have used that, as shown above, the maximum 
concurrence of rank 2 density matrices is its maximum 
eigenvalue. The upper bound in Eq. ([l3]) again coincides 
with C* for P4 = 0. 

In order to check whether C* is maximum even in gen- 
eral p4 cases, we have performed a numerical cal- 
culation whose scheme is similar to that in Ref. f|,||,U. 
We have generated 10,000 density matrices in a diagonal 
form with random four eigenvalues The maximum 
concurrence has been obtained among 1,000,000 density 
matrices generated by multiplying random unitary matri- 
ces in the circular unitary ensemble jl0| to each of 10,000 
diagonal matrices. The results are shown in Fig. |l| where 
the maximum concurrence are plotted as a function of 
the participation ratio (i?=l/Trp 2 ). 

When the density matrix is close to the pure state (R = 
1), the maximum concurrence is also close to unity, as 
expected. For i?>3, the states are always separable (Eq. 
(|l2|)) and the maximum is zero. In the region of 1 < R < 3, 
the maximum tends to decrease with the increase of R, 
but the points rather broadly distribute. The same data 
are plotted as a function of C* in the inset of Fig. |]. 
All points very closely distribute along the straight line 
of C = C* , and none of the points are present on the 
higher side of the line. This numerical result strongly 
support the hypothesis that C* gives an upper bound of 
the concurrence even in the general cases of P4^0. 

Accepting the hypothesis implies that all the states sat- 
isfying C*(pi)<0 become automatically separable. This 
condition of separability is looser than Eq. (fL2|). In fact, 
C*{pi) < is only a necessary condition of Trp 2 < 1/3. 
Difficulty of the rigorous proof of the hypothesis, if it 
is true, might relate to the difficulty of complete un- 
derstanding of separable-inseparable boundary in the 15- 
dimensional space of the density matrices due to its com- 
plex structure. We emphasize again that the numerical 
result strongly support that the hypothesis is true. 

It should be noted here that, when the eigenvalues of 
a density matrix satisfy a relation, C* is indeed maxi- 
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mum even for pi ^ 0. Any rank 4 density matrices are 
decomposed as 



p = (pi -P3 - 2vp2pI)|i>(i| + Pi, 



(14) 



where 1 1) is an eigenvector of p, and eigenvalues of p\ (not 
normalized) are {p3 + 2^/p2P4, P2, P3, Pi} ■ When eigenval- 
ues of p satisfies 



(15) 



P3 = Pi + Pi - \/PlPi, 



the purity of (normalized) p^ is equal to 1/3 and pi be- 
comes always separable. Therefore, using the convexity 
of the concurrence again, the upper bound of the concur- 
rence is proven to be C* for density matrices satisfying 
Eq. (|l5|). When P2=P3=Pi(< 1/4), M is reduced to the 
Werner state: 



M 



Pi] 1 ® 
1 - 



-><* 
Pi 



(|*+)(vl/ + | + |$-)($-| + |<I> + )($ + |), (16) 



whose eigenvalues satisfies Eq. (jlq). Therefore, EOF of 
the Werner states cannot be increased further by any 
unitary transformations. 

It is worth to test whether the states we propose have 
maximum entanglement in the other entanglement mea- 
sures. It has been shown that positive partial transpose 
is necessary condition for separability Q, and that it is 
also sufficient condition for 2x2 and 2x3 systems [ pd| . 
In 2 x 2 systems, when the density matrix is entangled, its 
partial transpose has only one negative eigenvalue [ jT2| . 
The modulus of the negative eigenvalue (En) is one of 
entanglement measures, and twice of En agrees with the 
negativity introduced in Ref. We have performed a 
numerical calculation similar to that for Fig. |l|, and ob- 
tained maximum of En for 10,000 random density matri- 
ces. In the inset of Fig. ||, those are plotted as a function 
of 

2E* N = max{0,2E N (j>i)} 



-P2 -Pi + VAPi -P3) 2 + (P2 - Pi) 2 , (17) 



which is the negativity of M. None of the points are 
present on the higher side of a straight line of En = E% 
as in the case of the concurrence. While it has been 
shown that two measures (EOF and negativity) do not 
induce the same ordering of density matrices with respect 
to the amount of entanglement [|| , the above numerical 
results suggest that M has the maximum amount of en- 
tanglement in both two measures. 

It should be noted here that E N = is equivalent to 
C* = 0, and therefore, the condition of C*(pi) < 0, for 
which all the states will be separable as mentioned be- 
fore, does not contradict to the condition of E N (pi) <0. 
Further, it will be natural to attribute the upper bound 
of entanglement, which is well described by C*(pi) and 
E N (pi), to the increase of the degree of the mixture of the 
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FIG. 2. The same as Fig. |lj but negativity is plotted. 

states. In this sense, each of C*(pi) and 2E* N (pi) (both 
distribute in the range of [—1/2,1]) can be considered 
as one of measures characterizing the degree of mixture 
such as the purity (or participation ratio), von Neumann 
entropy, Renyi entropy, and so on. 

Finally, as a simple application of the upper bound of 
EOF: 



E F (p) < H( 



1 + VI - C" 



). 



(18) 



we consider the situation generating entangled states by 
using the quantum gate consisting of two qubits, more 
concretely a controlled-NOT (CNOT) gate. In realistic 
situations, the coupling between the gate and its sur- 
rounding environment is inevitably present. The entire 
system consisting of the gate plus its environment hap- 
pen to be entangled by the coupling, and the mixture of 
the reduced density matrix describing the gate will be 
inevitably increased. 

In order to treat such decohered CNOT gates, we adopt 
the spin-boson model |l3|,[l4|], where each qubit is de- 
scribed as a spin i system, and the environment is ex- 
pressed as an ensemble of independent bosons. As the 
model of the CNOT gate, we choose the simplest Hamil- 
tonian: 



H G = - j(l - <Jcz) <8 <Jtx, 



(19) 



where c and t denotes the control-bit and target-bit, re- 
spectively. The state change after t = h/ (2R) corresponds 
to the change in the CNOT operation. In this paper, 
we demonstrate two types of gate-environment couplings. 
These are 

Hq E = ~u cz ^ B k{a\ + a k ), 

k 

H GE = -^(l - 0"cz) ® <?tx Bk( - Q l + Q k), ( 2 0) 
1 k 

where a k is an annihilation operator of a boson in the 
environment. Hq E describes the situation that only the 
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FIG. 3. EOF as a function of time in a decohered CNOT 
gate for several values of the coupling strength. uj c is the 
cut-off frequency of the Ohmic environmental mode, and j3 is 
the inverse temperature, (a) Hq E and (b) Hq E coupling. 

control-bit couples with the environment. Hq E may de- 
scribe the situation that the gate operation is achieved by 
irradiating a optical pulse which contains a noise coher- 
ent over the qubits as well as the pulse itself. For these 
phase-damping couplings, the time evolution of the re- 
duced density matrix describing the gate is analytically 
solved |l5| by assuming the product initial state for the 
entire density matrix: 



Ptot(0) = p(0) <g> p E , 



(21) 



where pe is the thermal equilibrium density matrix of 
the environment. Since we pay attention only to the 
generation of the entangled state, the initial state of the 
CNOT gate is chosen to be the pure state of (|0) c +|l) c )(8) 
|0) t - 

The time development of EOF for several values of 
the coupling strength K (oc ^2 k Bf,5{ui — LOk)) ar e shown 
in Fig. H for the (a) H^ E and (b) Hq E coupling. The 
values of the coupling strength K are chosen such that 
the values of the fidelity of the output state to the desired 
state in the absence of the decoherence are roughly 0.95, 
0.9 and 0.8, which are common to Figs. | (a) and (b). It 
is interesting to note that, while the fidelity is the same 
in two cases of coupling, the amount of the entanglement 
is significantly different. 

The upper bound of EOF (Eq. (|T8|)) for each coupling 
strength is shown by an arrow on the right side of each 
panel for comparison. Since EOF of a state has the phys- 
ical meaning of the asymptotic number of Bell pairs re- 
quired to prepare the state by using only local quantum 
operations and classical communication, comparing the 
difference in EOF will make sense. In Fig. | (a), the 
EOF is considerably lower than the corresponding upper 



bound, while EOF almost agrees with the upper bound 
in Fig. H (b). Therefore, with respect to the function gen- 
erating entangled states, the performance of the CNOT 
gate shown in Fig. || (b), is already optimal (or satu- 
rated) in a sense that there is no other way to increase 
further the amount of the entanglement than avoiding 
the increase of the mixture of the output density matrix, 
and thus, avoiding the decoherence itself. On the other 
hand, for the CNOT gate shown in Fig. |^ (a), there is 
room for improvement of the performance in principle, 
although we cannot show the detailed methods here. 

To conclude, we propose the mixed states in two 
qubits, which have a property that the amount of en- 
tanglement of these states cannot be increased further 
by any unitary operations. The property is proven when 
the rank of the states is less than 4, and when the states 
satisfy a special relation such as the Werner states. The 
results of the numerical calculations strongly support a 
hypothesis that these mixed states are indeed maximally 
entangled even in general cases. It will be extremely im- 
portant to verify the above hypothesis and to seek out 
the maximally entangled mixed states as well as the mea- 
sure in larger dimensional systems, for understanding the 
nature of entanglement of general mixed states and for 
the progress of the quantum information science and its 
applications. 
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